arXiv:1502.04816vl [astro-ph.EP] 17 Feb 2015 


Asteroid Models from Multiple Data Sources 


Josef Durech 

Astronomical Institute, Faculty of Mathematics and Physics, Charles University in Prague 


Benoit Carry 

Institut de Mecanique Celeste et de Calcul des Ephemerides 


Marco Delbo 

Lahoratoire Lagrange, UNS-CNRS-Observatoire de la Cote d’Azur 


Mikko Kaasalainen 

Tampere University of Technology 


Matti Viikinkoski 

Tampere University of Technology 

In the past decade, hundreds of asteroid shape models have been derived using the 
lightcurve inversion method. At the same time, a new framework of 3-D shape modeling 
based on the combined analysis of widely different data sources such as optical lightcurves, 
disk-resolved images, stellar occultation timings, mid-infrared thermal radiometry, optical 
interferometry, and radar delay-Doppler data, has been developed. This multi-data approach 
allows the determination of most of the physical and surface properties of asteroids in a single, 
coherent inversion, with spectacular results. We review the main results of asteroid lightcurve 
inversion and also recent advances in multi-data modeling. We show that models based on 
remote sensing data were confirmed by spacecraft encounters with asteroids, and we discuss 
how the multiplication of highly detailed 3-D models will help to refine our general knowledge 
of the asteroid population. The physical and surface properties of asteroids, i.e., their spin, 3-D 
shape, density, thermal inertia, surface roughness, are among the least known of all asteroid 
properties. Apart for the albedo and diameter, we have access to the whole picture for only a 
few hundreds of asteroids. These quantities are nevertheless very important to understand as 
they affect the non-gravitational Yarkovsky effect responsible for meteorite delivery to Earth, or 
the bulk composition and internal structure of asteroids. 


1. INTRODUCTION 


The determination of asteroid physical properties is an 
essential part of the complex process of revealing the nature 
of the asteroid population. In many cases, this process starts 
with obtaining observational data, continues with creating 
a model of the asteroid (i.e., its size, 3-D shape, and spin 
state, in the first approximation), and ends with interpreting 
new facts based on the model or a set of these. In this sense, 
modeling is a crucial mid-step between observations and 
theory. Results based on individual well-studied asteroids 
can be generalized to other members of the population. On 
the other hand, a statistically large sample of asteroids with 
known properties can reveal physical effects that play an 
important role for the whole population. 

In this chapter, we will build on the content of the Aster¬ 
oids III chapter by Kaasalainen et al. ( 2002ah about aster¬ 
oid models reconstructed from disk-integrated photometry. 
Although visual photometry still remains the most impor¬ 


tant data source of the modeling, the main progress in this 
field since Asteroids III has been the addition of comple¬ 
mentary data sources. Many of these data sources are disk- 
resolved, thus containing much more information than disk- 
integrated data. This shift in paradigm - using photome¬ 
try not alone but simultaneously with complementary data 
- was mentioned in the last paragraph of the Asteroids III 
chapter as ‘'perhaps the most interesting future prospect”, 
and we are now at this stage. In the following, we will 
review all data types suitable for inversion, their sources, 
uncertainties, and how they can be used in modeling. 

When describing the methods of data inversion and the 
results obtained by these methods, it is also important to 
emphasize caveats, ambiguities, and possible sources of er¬ 
rors. Although the description of what can be obtained from 
different data sources is exciting, the knowledge of what 
cannot, i.e., what are limitations of our data sets, is of the 
same importance. Omitting this may lead to over interpre¬ 
tation of results. 
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This chapter is structured as follows. First we review 
the main principles of the multimodal inverse problem in 
Sect. |2] Then, in Sect. [3 we discuss each data type and 
their contribution to model characteristics and details, and 
we describe some extensions of the predominant model. In 
Sect, m we discuss the main results based on lightcurve in¬ 
version and multimodal asteroid reconstruction. We con¬ 
clude with prospects for the future in Sect.|5] 


2. THEORETICAL ASPECTS OF INVERSION AND 
DATA FUSION 

Asteroid physical model reconstruction from multi¬ 
modal data is, by its nature, a mathematical inverse prob¬ 
lem. It is ill-posed; i.e., the uniqueness and stability prop¬ 
erties of the solution are usually not very good unless the 
data are supported by a number of prior constraints. Fur¬ 
thermore, it is not sufficient just to fit some model to the 
data numerically and try to probe the solution space with 
some scheme. Although there are more approaches to the 
problem of asteroid shape reconstruction, they are usually 
dealing with only one data type and we mention them in the 
next section. Here, we describe the problem in a general 
way in the framework of generalized projections: our data 
are various 1-D or 2-D projection types of a 3-D model, and 
understanding the fundamental mathematical properties of 
the inverse projection mapping is essential. This includes 
a number of theorems o n uniqueness, information content , 
and stability properties iKaasa lainen and Lambe t^ 2006t 
\Kaasalainen , 2011 : Viikinkoski and Kaasalainen . 20141) . 


Let the projection point in the image plane (plane- 
of-sky or range-Doppler) of the point Xq on the body be 
mapped by the matrix A: Define the set I(^) 

for any ^ as 

AO = {x|5(C, x; R, t) h{x- M, R, t) = 1}, (1) 

where we have explicitly shown the time t and the ad¬ 
justable parameters: M for the shape and R for the rotation. 
The projection point function x) = 1 if A{R, f)x = 
and p = 0 otherwise. The ray-tracing function h = 1 if x is 
visible (for occultation, thermal, and radar data), or visible 
and illuminated (for disk-resolved imaging and photometry 
in the optical); otherwise h = 0. The set 1(0 is numerable 
and finite. The number of elements in I(^) is at most one 
for plane-of-sky projections (each point on the projection 
corresponds to at most one point of the asteroid’s surface); 
for range-Doppler, it can be more (more points on asteroid’s 
surface can have the same distance to the observer and the 
same relative radial velocity). Generalized projections, i.e., 
all the data modes presented in Sect. [3 can now be pre¬ 
sented as scalar values p(^) in the image field H: 

piOt) = j fiOv) E S{x; M, R, L,t) dr], (2) 

xei(77) 

where L denotes the luminosity parameters (for scattering 
or thermal properties), and the luminosity function is de¬ 
noted by S. The function / is the point-spread, pixellation. 


or other transfer function of the image field. For interfer¬ 
ometry, it is typically the Fourier transform kernel. In fact, 
the reconstruction process works efficiently by taking the 
Fourier transfor m of any image type rather than usi ng the 
original pixels iviikinkoski and Kaasalainal. 12014 ). For 
lightcurves, / = 1 (and ^ is irrelevant, p(^) is constant). 
The surface albedo is usually assumed to be constant, al¬ 
though its variegation can be included in S by the param¬ 
eters L if there are high-quality disk-resolved data. In the 
case of lightcurves only, we can get an indication of non- 
uniform alb edo and compensate for th is with a (non-unique) 
spot model iKaasalainen et al . , 200 ih . 

The multimodal inverse problem can be expressed as fol¬ 
lows. Let us choose as goodness-of-fit measures some func¬ 
tions 6i, i = 1,..., n, of n data modalities. Typically, <5 is 
the usual form between Pmodei and Pobs- Our task is 
to construct a joint 6tot with weighting for each data mode: 


Stot{P,D)=Si{P,Di)+J2x,_iS,{P,D,), D = {D,}, 

( 3 ) 

where Di denotes the data from the source i, Ai_i is the 
weight of the source i, and P = {M, R, L} is the set 
of model parameter values. The best-fit result is obtained 
by minimizing (5tot with nonlinear techniques, typically 
Levenberg-Marquardt for efficient convergence. Regular¬ 
ization functions r{P) can be added to the sum; these con¬ 
strain, for instance, the smoothness of the surface to sup¬ 
press large variations at small scales, the deviation from 
principal-axis rotation to force the model to rotate around 
the shortest inertia axis (as suming uniform density), or 
the g ravitational slope, etc. iKaasalainen and Viikinkoski . 
20i3). 


The modality (and regularization) weights Ai are deter¬ 
mined us ing the maximum compatibility estimate (MCE! 
princ iple iKaasalaineA 201 it Kaasalainen and Viikinkoski , 
20121) . This yields well-defined unique values that are, in 
essence, the best compromise between the different data 
sets that often tend to draw the solution in different direc¬ 
tions. Moreover, MCE values of weighting parameters are 
objective, not dependent on users choice, although their 
values are usually close to those determined subjectively 
based on experience. Plotting various choices of weights 
typically results in an L-shaped curve shown in Fig. [3 the 
best solution is at the corner of the curve. In this way, the 
reconstruction from complementary data sources is possi¬ 
ble even if no single data mode is sufficient for modeling 
alone. 

For practical computations, the surface is rendered 
as a polyhedron, and S and h are computed accord¬ 
ingly with ray-tracing iKaasalainen et ali 120011) . Rather 


than using each vertex as a free parameter, the surface 
can be represented in a more compact form with spheri¬ 
cal harmonics series (for starlike or octantoid; i.e., gen- 
eralized starlike, shapes) or subdivision control points 


^Kaasalainen and Viikinkoskn.Wdl'MViikinkoski et a/.Ll2015h 
These shape supports are essential for convergence: they 
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Fig. 1.— The level of fit for lightcurves and occultation data 
for different weighting between the two data types. The optimum 
weight is around the “comer” of the L-curve. Each dot corre¬ 
sponds to an inversion with a different weight A. 



Fig. 2.— Model of (41) Daphne reconstmcted from lightcurves 
and adaptive-optics images using subdivision surfaces (top) and 
octantoids (bottom). The general shape remains stable, even if 
small-scale features slightly change. 


allow flexible modifications of the surface with a moderate 
number of parameters while not getting stuck in local min¬ 
ima or over-emphasizing the role of regularization functions 
when searching for the best-fit solution. Each shape sup¬ 
port has its own characteristic way of representing global 
and local features. For example, the octantoid parametriza- 
tion 


r x{0,ip) = sin0cost/?, 

x(0, (p) = } y{e, p) = sin B sin </?, 

\ zle,ip)= 

(4) 

where a, b, and c are linear combinations of the (real) spher¬ 
ical harmonic functions p), with coefficients {aim}, 

[bim], and [cim], respectively, is easy to regularize glob¬ 
ally while retaining the ability to produce local details. The 
coordinates {9, </?), O<0<7r, 0<(/3< 2tt, parametrize 
the surface on the unit sphere S'^ but do not represent any 
physical directions such as polar coordinates. 

This inverse problem is a typical example of a case 
where model and systematic errors dominate over random 
measurement errors. Thus the stability and error estima¬ 
tion of the solution are best examined by using different 
model types (Fig. |2]i. In the case of shape, for example, 
the reliability of the features on the solution can be checked 
by comparing the results obtained with t wo or more shape 


suppo rts (starlike, octantoid, suhdiyision: \Viikinkoski et al. 


2015h . This yields better estimates than, e.g., Markov chain 


Monte Carlo sequences that only investigate random error 
effects within a single model type. 

A particular feature of the model reconstruction from 
disk-resolved data is that the result is dominated by the 
target image boundaries rather than the pixel brightness 
distribution within the target image. This is because 
the information is contained in the pixel contrast which 
is the largest on the boundary (occultations are special 


cases of this as they are samples of the boundary con¬ 
tour). This is very advantageous when considering the 
effect of model errors in luminosity properties (scattering 
or thermal models); it is sufficient to have a reasonable 
model, and the result is not sensitive to the parameters 
L. Thus, for example, Atacama Large Millimeter Array 
data can be used for efficient reconstruction even with 
a very approximate semianalytical Fourier-series thermal 
model - more detailed mod els have hardly any effect on 
the shape solution {S&ct. [i .5\\Viikinkoski andKaasalainen . 


2014 : Viikinkoski et al. . 2015b . 


3. DATA AND MODELING 

We describe all data types that can be used, the way of 
collecting the data, their accuracy, typical number of aster¬ 
oids for which the data exist, and expectations for the fu¬ 
ture. We also discuss a typical result of inversion - what 
is the resolution of the model and how many targets can be 
modeled (Table [Til. 


3.1. Photometry 

Disk-integrated photometry is, and will always be, the 
most abundant source of data, because it is available for 
essentially every single known asteroid. Because aster¬ 
oid brightness periodically changes with its rotation, fre¬ 
quency analysis of asteroid lightcurves provides asteroid ro¬ 
tation periods - the basic physical property derivable from 
time-resolved photometr y. The regular l y upd ated Asteroid 
Lightcurve Database of \Warner et al\ (l2009b available at 
http://www.minorplanet.info/lightcurvedatabase.html now 
contains rotation periods and other physical parameters for 
almost 7000 objects, for about half of them their rotation 
period is secure and unambiguous. The role of amateur as¬ 
tronomers in this field is traditionally strong, getting even 
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stronger with increasing level of their technical and soft¬ 
ware equipment. Hundreds of asteroid lightcurves are pub¬ 
lished quarterly in the Minor Planet Bulletin; most of them 
are then archived in the Asteroid Light Curve Database 
at the Minor Planet Center site (http://mpc.cfa.harvard. 
edu// light_curve2/l i eht_cu rve.php) in the ALCDEF stan¬ 
dard iWarner et al\. 1201 ih . The efficiency of lightcurve 


production can be increased by dedicated wide-held photo- 
metric surveys ( Masiero et al . , 20091: Polishook et al . , 20121 
for example), although the period determina tion from un 


dersa mpled lightcurves is often ambiguous iHarris et al. 

20121 ) . 


For period determination, a single lightcurve cover¬ 
ing the full rotation is sufficient. However, a set of such 
lightcurves observed at different geometries (asteroid il¬ 
luminated and seen from various directions) is needed to 


lightcurve inversion method of 

Kaasalainen and Torppa 

(2001): Kaasalainen et al. (2001 

) was already reviewed 

in Asteroids III ( 

Kaasalainen et al. 2002al). Since then. 


the method has been widely used and hundreds of aster¬ 
oid models have been derived. They are publicly avail¬ 
able at the Database of Asteroid Models from Inver¬ 
sion Techniques (DAMIT, http://astro.troja.mff.cuni.cz/ 
projects/asteroidsSD). The reliability of the method was 
proved by comparing its results with independent d a ta suc h 
as laboratory asteroid model {Kaasalainen et ail 120051) . 


adantive-o ptics images iMarchis et al.l 120061) . stellar oc- 


cultations iPurech et al.i \2Qllh. or space craft images of 


asteroids (2867) Steins iKeller et al.[ 2010l) and (433) Fros 
( Kaasalainen et al. . 2002ah . 

From disk-integrated photometry alone, only a global 
shape without any small-scale details can be derived. Be¬ 
cause the reflectivity of the surface is not known, the 
models are not scaled and the information about the size 
has to come from complementary data. To avoid over¬ 
interpretation and artifacts of the modeling, the shapes 
are usually represented by a convex model. This allows 
to work not in the obvious radius parameter space but 
in the Gaussian image space (describing a convex body 
by the curvature of its surface). This is less intuitive 
but it makes the inverse problem less vulnerable to er- 
rors of data and model because of the Minkowski stability 
[ Lambers and Kaasalaineri 2001 ) - even if the areas of in¬ 
dividual surface facets may change significantly for slightly 
different data sets, the global convex shape changes very 
little. 

From the practical point of view, we are interested in 
finding a unique solution of the inverse problem. To guaran¬ 
tee this, observations covering a sufficiently wide range of 
viewing and illumination geometries are needed. For a typi¬ 
cal main-belt asteroid, it means observations during several 
apparitions. For a close-approaching near-Farth asteroid, 
several months could be sufficient. For more distant objects, 
we can in principle observe changing lightcurve amplitude 
because of changing aspect, this would anyway take many 
decades. Nevertheless, trans-Neptunian objects (TNOs) 



Fig. 3.— Sparse-in-time photometry of (21) Lutetia obtained at 
US Naval Observatory plotted as a function of time (top) and phase 
angle (middle). The brightness was reduced to the unit distance 
from Earth and Sun. The scatter in the phase plot is caused by 
the intrinsic noise of the data and the rotational and aspect effects. 
The bottom panel shows residuals for the best-fit model (Fig. 
plotted over the phase curve (solid curve). 


and Centaurs can be never observed from Farth at geometry 
significantly different from opposition, wh ich is not suffi - 
cient to reconstruct a unique convex model ( RusselL 1906h . 
If the lightcurves observed in various filters are different, 
there is some spectral/color variegation ov er the surface. 


and a crude color map can be reconstructed iNathues et al. 


20Q.5h. 


A typical outcome of lightcurve inversion is a convex 
shape model that describes the global characteristics of the 
real asteroid. It is also easy to obtain nonconvex versions 
with the general inversion procedures (see Sect. El), but then 
one should produce several solutions with various model 
types and parameters, and be very cautious about the re¬ 
sults ( Viikinkoski et al\. 2015 ). In general, disk-integrated 
photometry contains very little information about noncon¬ 
vexities unless they are very pronounced or observed at 
very high phas e angles where shadowing effe cts play an 
important role jPurech and Kaasalaineri l2003l) . Noncon¬ 
vex models seldom fit lightcurves better than convex ones 
simply because the latter typically already fit the data down 
to noise level. This sets the resolution limit of photome¬ 
try (see the discussion in Kaasalainen et al[ 2001 , 2002a ; 


Viikinkoski et a/.Ll2015L and references therein). 


The relative accuracy of the sidereal rotation period de- 
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to the noise level [DurechetaL, 2005h . However, the total 


Fig. a The best-fit model based on this data set has the 
same rotation period and a similar spin axis direction as the 
models based on much larger and multiple data sets (Fig. ISJ. 

When using noisy sparse data or only a limited number 
of data points, a simple shape model of a triaxial ellipsoid 
is usually sufficient to model the data and to derive the cor¬ 
rect period and spin vector orientation ( Cellirwet ^ , 20091 : 


10 15 20 

Period [hours] 

Fig. 4.— Period search for (21) Lutetia. Each point represents 
a local minimum in the parameter space to which the optimiza¬ 
tion algorithm converged. The lowest RMS residual (arrow) cor¬ 
responds to the best-fit model shown in Fig.|^ 


termined from lightcurves is of the order of 10“^ or better, 
depending mainly on the time span of observations. The 
direction of the spin axis can be determined with an ac¬ 
curacy from a couple of degrees for the models based on 
many decades of observations with many dense lightcurves, 
to more than twenty degrees for models based on limited 
and noisy data. For asteroids orbiting in the ecliptic plane, 
the geometry for an Earth-based observer is limited to that 
plane. Then the disk-integrated brightness of a body with 
surface described by radius vector (x, y, z) and pole di¬ 
rection (A, /?) in ecliptic coordinates is the same as for a 
body (x,y,—z) with the sp i n axis direction (A -I- 180°,/3) 
iKaasalainen and Lambem 120061) . That is why for a typ¬ 
ical main-belt asteroid, there are usually two equally good 
mirror shape solutions with about the same pole latitude and 
pole longitude difference of about 180°. This ambiguity 
can be removed with disk-resolved plane-of-sky projections 
(e.g., images, stellar occultations, see Sect. 13.2113.41 i. 

Apart from “classical” lightcurves where the sampling 
of brightness variations is dense with respect to the rota¬ 
tion period, there are also data that are sparse in time. Such 
data sets typically consist of only a single to a few bright¬ 
ness measurements per night. Provided the whole sparse 
data set is internally calibrated, it can be used the same 
way as a standard lightcur ve that would be ext remely long 
and very sparsely sampled ( Kaasalainen . 20o3) . These data 
are routinely provided by all-sky astrometric surveys with 
various - usually poor - photometric quality. Current sur¬ 
veys provide data of accuracy scarcely better than 0.1 mag. 
Given a typical lightcurve amplitude of 0.3 mag, the sig¬ 
nal is often drowned in noise and systematic errors. This 
leads to many phy sically acceptable m odels fitting the data 


Cellino and Dell’Ord. 2012t Carhoenani et al. . 2012 ). The 

adyantage of this approach is that the shape is described 
by only two parameters (axes ratios) and the scanning of 
the period parameter space is much faster than with general 
shapes. 

The on-going (Pan-STARRS, Catalina, Gaia) and future 
(ATLAS, LSST) suryeys will proyide new data eyery night 
for essentially eyery known asteroid. Using this data for au¬ 
tomated lightcurye inyersion with well-mastered treatment 
of systematic effects, recognition of the best-fit models, def¬ 
inition of uniqueness of the solution etc., is the main chal¬ 
lenge for the future lightcurye inyersion. 

3.2. Remote sensing disk-resolved images 

The most direct way to obtain information on the shape 
of an asteroid is to take pictures of it. The apparent shape as 
yisible on the plane of the sky is delimited by the limb and 
the terminator, and multiple yiews obtained while the tar¬ 
get rotates can fully characterize its 3-D shape. To resolye 
the small angular diameter (< 0.5") sustained by asteroids, 
large facilities are, howeyer, required. 

In the 1980s, speckle imaging or speckle interferome¬ 
try proyided “f/ze first glimpses of an as teroid’s surface” 


(i.e., (4) Vesta by Drummond et al. . 1988h . This technique 


is based on the analysis of the speckle pattern in the images 
of astronomical sources obtained through large telescopes 
at high magnification power and yery short integration time. 
The aim is to oyercome the blurring effects due to the astro¬ 
nomical seeing and to attain diffraction-limited resolution 
images. Speckle interferometry has been commonly used 
to study t he size, shapes, and surface features of the largest 


asteroids (^Drummond et ali I1985[ 19881: Raeazzoni et al. . 
2 OOOI: Cellino et al\. 20031) 7 With the launch of the Hub- 


amount of data is huge and it has been shown that at least 
for some asteroids, models can be successfully deriyed from 
these data or from their combination wit h dense lightcuryes 
i Purech et ali 20091 : Hanus et al. . 201 ih . 

As an example, we show sparse photometry for asteroid 
(21) Lutetia from US Nayal Obseryatory (260 points coy- 
ering ten years) in Fig. [3 and the period search results in 


ble Space Telescope (HST) in orbit and the first light of the 
large (10-m class) ground-based telescopes equipped with 
adaptiye-optics fed cameras (e.g., W. M. Keck, European 
Southern Obseryatory (ESO) Very Large Telescope (VLT), 
Geminis, Subaru), the importance of speckle interferometry 
has decreased. 

The critical issue in direct imaging is of course the an¬ 
gular resolution. Any image is the result of the conyolu- 
tion of the object on the plane of the sky with the instru¬ 
ment response, the point spread function (PSE). In space, 
the PSE is stable and corresponds to the diffraction pat¬ 
tern of the telescope. Erom the ground, the atmospheric 
turbulence constantly deforms the PSE and blurs the im¬ 
ages, hence the need of real-time correction of the PSE 
by adaptiye optics (AO). The technical challenges of send¬ 
ing a large telescope in space and of building deformable 
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Fig. 5.— Deconvoluted adaptive optics image (left) of asteroid 
(41) Daphne, and the corresponding image of the reconstructed 
model (right, projection of the the octantoid model from Fig. IS. 
The model is shown under artificial illumination that enhances 
its 3-D shape and is different from the simple light-scattering 
model used for the inversion. The AO image has resolution of 
lOmas/pixel. The model was reconstructed from 14 AO images 
obtained with ESO VLT and several lightcurves. 


mirrors, explain why the first disk-resolved images in the 
1990s were still limited in resolution and only the largest 
asteroids (1) Ceres, (2) Pallas, and (4) Vesta have been im¬ 
aged. In the decade since Asteroids III, numerous stud- 
ies have been based on direct imaging of asteroids (e.g. , 
Thomas^etal^T^S, 2008; Descamos et al. 


2008t Schmidt et a/.l 20091: Marchis et al . 2013 ). 

Both from space or with adaptive optics on the ground, 
however, the contrast and angular resolution can be im¬ 
proved by deconvolution of the image by the PSF. This 
is particularly true for images acquired from the ground 
with residuals from non-perfect AO correction. Decon¬ 
volution is an ill-posed problem, but robu st algorithms 


adapt e d to planetary image s are available iConan et al.^ 
2 OO 0 I: Musnier et ai . 20041: Homet aL . 200?!) and have 


been validated on sky ( Witasse et al. . 20061) . It is never¬ 
theless the most critical part of the post-processing, as an 
incorrect deconvolution can introduce a systematic error on 
the apparent size. An example of an AO image and the re¬ 
constructed shape model for asteroid (41) Daphne is shown 
in Fig.|5] 

Current facilities deliver an angular resolution of about 
30 to 50 milliarcseconds (mas), depending on the wave¬ 
length. The apparent shape can therefore be measured for 
asteroids with an apparent diameter larger than « 80-100 
mas, i.e., a couple of hundred targets. Simulations and ob- 
servations of known targets such as the satellites of Saturn 


{ Marchis etal. . 20061: Drummond et al. . 2009; CarrX 2009) 
have shown that a precision of a few mas can be derived on 
the 2-D profile on the plane of the sky, corresponding to 
only a few kilometers for main-belt asteroids. With upcom¬ 
ing large telescopes (30-)- m, such as the Thirty Meter Tele¬ 
scope or the European Extremely Large Telescope ), the an¬ 
gular resolution will be improved by a factor 3^, providing 
more than 500 targets. Second generation instruments with 
extreme AO foreseen on these telescopes should allow the 
observation of about 7000 as t eroids , with sizes of only a 
few kilometers { Merline et g/.L 2013 ). 



Eig. 6.— Different shape models for (21) Lutetia derived 
from different sets of data. The level of detail increases with 
more data from to p to bottom: s parse photometry (Sect. [LB . 
dense lightcurves iTorvoa et g/.L 12003 ), lightcurves and AO 
j Carry etg/.L 12010^ . Rosetta flyby ( Sierks et ^1201 ih . 


Using disk-resolved i mage s , giant craters have b een dis¬ 


covered ((Thomas et alx Il997t \Conrad et al. 


guity in spin solu t ions ha ve been solved (Marchis et al. 


2OO7I). ambi- 


20061 : Carry etal. . I20 1 0a), albedo m a ps hav e been con¬ 
structed { Li et al. . 12006 : Carry et al. . 2008 ). convex 3- 
D shape models from lightcurves have been set to scale 
( Hanus et al ^ 2013al). and full 3-D shape mode ls deter¬ 
mined ( Carry et al. . 201 Obi: Descamos et al. . 2011 ). 


3.3. Flybys 

Asteroid imaging and mapping in situ is essentially an 
extreme example of disk-resolved images. The model¬ 
ing in such cases is more a cartographic than an inverse 
problem as the data are abundant, directly usable (contain¬ 
ing identification points on the surface, etc.), and high- 
resolution (there is no ill-posedness). However, asteroid 
flybys, during which roughly a half of the target is likely 
to remain unseen (not visible and illuminated), pose a spe¬ 
cial problem: how to see the dark side? The principle 
here is to use the high-resolution map of the seen side. 
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constructed with a number of metho ds of photoclinometry , 
photogrammetry, and image fitting ( ^eusker et al\ . 12012 ; 


Gaskell et ail l2008t Uorda et alV 2012h . as a constraint in 
the otherwise same multi-mode inverse problem as with 
gr ound-based observations. The pro cedure is described 
m lKaasalainen and Viikinkoskn (120121) : with it, the recon- 
struction of the da rk side, such as those of (2867) S teins 
{ Keller et al ]. l2010h or (211 Lutetia \Sierks et aZ.I . Eoi l[ see 
also Fig. is more detailed than from, e.g., lightcurves 
alone. This is because half of the target is accurately re¬ 
constructed, with practically no error margin: therefore the 
fluctuation margin of the dark side is considerably smaller 
as well. 


3.4. Stellar occultations 

The observation of a stellar occultation consists in 
recording the duration of the disappearance of a star behind 
the asteroid. Knowing the apparent motion of the asteroid 
on the plane of the sky, obtained from its ephemeris, this 
duration can be converted in a physical length on the disk 
of the asteroid, called a chord. Provided several observers 
record the same event from different locations on Earth, 
the 2-D profile of the asteroid i s drawn on the plane of the 
sky { Minis and Dunham . 19891) . The main difference with 
disk-resolved imaging resides in the profile, made by only 
the limb for occultations, and limb-tterminator for imag¬ 
ing. In fine, both techniques provide the 2-D profile of the 
target as projected on the plane of the sky at the epoch of 
observations (Fig.|7]i. 

Disk-resolved imaging and stellar occultation are, how¬ 
ever, radically different in term of facilities, data process¬ 
ing, potential targets, reproducibility, and achievable preci¬ 
sion. For stellar occultations, the properties of the occulted 
star matter generally more than the actual target: the aster¬ 
oid. If the occulted star is bright enough, its occultation, 
even by a very small asteroid, can be detected with small 
aperture telescopes. This is of course the main advantage 
of stellar occultations, where the apparent size and shape of 
potentially any asteroid can be measured. Moreover, this 
technique can be successfully used also for distant TNOs, 
whi ch have angular sizes too small to be resolved by imag¬ 
ing { Sicardv et al . 201 3 ). In practice, however, a given as¬ 
teroid will only seldom occult bright stars. Measurements 
are thus hard, if not impossible, to reproduce. 

Stellar occultations are nevertheless extremely valuable. 
The accuracy of the timing is dictated by time-series pho¬ 
tometry, and can therefore be extremely precise. An un¬ 
certainty of 50 ms in timings converts into only a 300 m 
uncertainty in the length of the chord, typical for a main- 
belt asteroid (at 1.5 au from Earth with an apparent motion 
of lO'Vh). The main source of uncertainty is, however, the 
absolute timing of each chord, required to align them on 
the plane of the sky. Most of historical occultations were 
recorded by naked eye, and suffered from this. Since a 
decade ago, thanks to the availability of low-cost position¬ 
ing and timing systems (e.g., GPS), stellar occultations are 



Fig. 7.— An example of a well-covered occultation that pro- 
vides almost an “imag e” of the projected asteroid (135) Hertha 
iTimerson et al\. 120090 . The model of Hertha based on its 
lightcurves and occultation data is over plotted. The dashed line 
is a negative observation constraining the northern part of the 
model’s projection. Formal timing uncertainties are between 0.02 
and 0.2 s, corresponding to 0.12-1.2 km in the projection, which 
is below the resolution of the model. The RMS of the fit is 1.9 km. 


being more and more valuable. 

The main contributors to the field are currently amateur 
astronomers: for a given event, observers have to move to 
set themselves on the predicted occultation path on Earth. 
Small aperture (< 20 cm) mobile stations are therefore ideal 
for recording stellar occultations. Because of the uncertain¬ 
ties of the star and asteroid positions on sky, there is gener¬ 
ally an uncertainty of a few tens of kilometers in the loca¬ 
tion of the occultation path on Earth, requiring observers to 
spread over large area to cover the event. This usually pro¬ 
hibits dense coverage of the asteroid profile. When only a 
couple of chords are available, the event provides only lim¬ 
ited information on the size, if any on the shape. Current 
occultation predictions of sufficient accuracy concern only 
stars in the Hipparcos catalogue and large (at least tens of 
kilometers) asteroids. With the upcoming publication of the 
ESA Gaia stellar catalog, and update of asteroid orbits, this 
position uncertainty is expected to drop signihcantly and 
future occultations will b e easier to predict, hence observe 
( Tanea and Delho . 2007h . 

Overall, stellar occultations can provide precise mea¬ 
surements of the size and shape of an asteroid, as projected 
on the plane of the sky. However, events are rare for a given 
target. Occultations are therefore very valuable in com¬ 
bined data sets, as for instance, to set scale to otherwise 


dimensionle s s 3-D shape models (&.s... \Durech et a/.Ll201 1 


Harms et alx. 1201 Ih. From almost 2500 occultations com- 
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piled hv \Dunham et al\ (120141) . there are about 160 “good” 
ones that allow a reliable determination of asteroid’s size 
and about 40 “excellent” ones that show details in asteroid’s 
profile. 

3.5. Interferometry 

Another technique to overcome the limitations of small 
angular sizes of asteroids, in order to measure their sizes, 
shapes, and possible presence of satellites, is interferome¬ 
try. An astronomical interferometer combines coherently 
(i.e., conserving the phase information) the light from two 
or more apertures of the same telescope or of distinct tele¬ 
scopes spaced by a distance B. The spatial resolution (in 
radians) is of the order of X/B where A is the wavelength. 

In the following, we give a basic introduction to interfer¬ 
ometry of aster oids. Further d etails can be found in the fol¬ 
lowing works: Janko^\ ( 2010l) and Matter et al. ( 2013 ) for 
a broad introduction to astron omical optical interf erometry, 
methods and instrumentation; Delho et al. ( 2009h for a de¬ 
scription of the techniques, and the models for deriving the 
size and basic shape proprieties of asteroids from the ESO 
Very Large Telescope Interferometer ( VLTI) MID-infrared 
Interferometric instrument (MIDI) data; CatT^_et^ (12015 ) 
for an extension of the technique o uDelho et al\ (120091) to 
the deter mination of the sizes and the separation of binary 
asteroids; Matter et al\ ( 201 ll) for a description and the ex¬ 
tension of a thermophysical model to the analysis of inter¬ 
ferometric data of asteroids with the aim of obtaining sur¬ 
face properties such as the thermal inertia. 

Interferometers measure the coherence function of the 
source, also called the interferometric visibility, which is 
given by the ratio between the correlated and the total 
flux. The correlated flux is the amount of flux in the in¬ 
terferometric fringes. More precisely the correlated flux is 
the Fourier transform of the brightness distribution of the 
source measured on the int erferometer’s basel ine projected 
on the plane of the sky (see Delho et 'ai\. 20091 for example, 
and references therein). 

The Fine Guidance Sensors (FGS) aboard the Hubble 
Space Telesco pe (HST) are optica l white-light shearing in¬ 
terferometers iNelan et all 119981) that combine the light 
from distinct apertures of the HST primary mirror and have 
been used to measure the size and basic shape proper - 
ties of asteroids iHestroffer et a/.Ll2002l : lranga et adl2003l) . 
HST/FGS data have clearly demonstrated the bi-lobed na¬ 


ture o f some bodies such as (216) Kleopatra iTanea et al. 


20011) and (624) Hektor iTansa et ail l2003h . but were not 
able to de tect the presence of the little moons orbiting these 
asteroids iPescamos et al.Ll201 luMarchis et a/.Ll2014l) . due 
to the large magnitude difference between the asteroids and 
their respective moons. Despite the impressive angular res¬ 
olution of a few mas, corresponding to a few kilometers 
at a distance of 1.2-1.5au, a clear limitation of HST/FGS 
asteroid studies is the bright limiting magnitude of the in - 
strument of about V ^ 13-14 mag uanea et alX. l2003h . 
A recapitulation of the HST/FGS asteroids results can be 



Angular Diameter (mas) 



Diameter (km) 

Fig. 8.— Visibility and total flux for an asteroid at 2.5 au 
from Sun and 1.5 au from Earth as a function of its diameter 
for different values of the baseline. The MIDI limiting cor¬ 
related flux (flux times visibility) is 0.5-1 Jy and 10-20 Jy 
for the UTs and the ATs respectively. Similar limiting fluxes 
are expected for MATISSE. See text for further details. 


found x mDotto et al\ (120021) . 

Ground-based interferometry is limited to observation 
of bright targets because of the requirement to take ex¬ 
posures shorter than the atmospheric coherence time of a 
few milliseconds. The limiting magnitudes of ground-based 
long-baseline interferometers and the intrinsic low surface 
brightness of asteroids, have prevented the use of these in¬ 
struments for the studies of small solar system bodies until 
the availability of the VLTI. 

Interferometry in the mid-infrared was proven to be also 
very sensitive to the global shape of asteroids and als o to 
their surface characteristics { Matter et g/.L 201 ll 2013 ). In 
particular, since observations are typically carried out in the 
thermal infrared (8-13 pm), MIDI data are sensitive to the 
surface temperature distribution, which is strongly affected 
by the value of the thermal inertia {Delho et al, this vol¬ 
ume). 

Figure[8]shows the visibility and the total flux as a func¬ 
tion of the angular diameter of an asteroid at 2.5 au from 
Sun and at 1.5 au from Earth, where 1 mas roughly corre- 
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spends to 1 km on the asteroid. Note that main-belt aster¬ 
oids smaller than ~20km cannot be easily observed with 
MIDI at the VLTI. There are about thousand known aster¬ 
oids with diameter above 20 km, implying that interferom¬ 
etry is potentially an interesting technique for shape mod¬ 
elling. On the other hand, as interferometry in the thermal 
infrared is sensitive to the spatial distribution of the tem¬ 
perature on the asteroid surface, this technique can be used 
to determine thermophysical properties of as teroids when 


the body sha pes are known as demonstrated hv xMatter et al 
(I201lil2013b . 

Each VLTI baseline can be used with MIDI only one at 
the time. MIDI will likely be decommissioned in the near 
future to be substituted after 2016 by the Multi AperTure 
mid-infrared SpectroScopic Experiment (MATISSE). This 
instrument will combine up to four Unit Telescopes (UTs) 
or Auxiliary Telescopes (ATs), allowing six simultaneous 
baselines. This feature will enable us to measure the spa¬ 
tial distribution of the infrared flux along different direc¬ 
tions. MATISSE will also measure closure phase relations 
and thus offer an efficient capability for image reconstruc¬ 
tion. In addition to the N band, the MATISSE will also 
operate in the L and M bands. Unfortunately, MATISSE is 
not expected to be more sensitive than MIDI. 

Another second generation instrument at the VLTI, 
GRAVITY, that will combine the light from all four UTs, 
will offer further improvements in spatial resolution com¬ 
pared to MIDI and MATISSE (though with more a more 
stringent V < 11 limiting magnitude). It will provide near- 
infrared adaptive-optics assisted precision narrow-angle 
(about 4") astrometry at the 10 pas level in the K band 
(2.2 pm). Both the reference star and the science object 
have to lie within the ~ 4" field of view. In imaging mode, 
GR AVITY can achieve a res olution of ~ 3 mas in the near- 
IR iEisenhauer et al\. l2008h . The imaging mode can be 
interesting to precisely measure the sizes and the orbits of 
the satellites of large asteroids, the latter with V < 11 mag. 

AMBER is the current near-infrared focal instrument of 
the VLTI. It operates in the J, H, and, K bands (i.e., from 
1.0 to 2.4 pm). The AMBER limiting magnitude for aster¬ 
oid observations is U ~ 9 mag. There is only a handful of 
asteroids brighter than this limit. These bodies are also the 
largest ones with angular extensions generally >100 mas, 
implying a very low visibility in the J-H-K with the UTs 
baselines. Although the photometric flux of these few aster¬ 
oids is such that If < 7.5 mag, their correlated magnitudes 
due to the low visibilities are much higher than 7.5 mag, 
preventing their fringe detection and tracking in the near- 
IR. 

However, one of the most interesting instruments for as¬ 
teroid ground based interferometry is the Large Binocular 
Telescope Interferometer (LBTI). It consists of two 8.4-m 
telescopes mounted side by side in a single mount, with 
a 14.4-m center-to-center spacing. This configuration of¬ 
fers a unique capability for interferometry of a Eizeau beam 
combination. This offers a wide field of view (^10-20") 
and low thermal background. Eor example, the LING- 



arcsec' 


Eig. 9.— Simulated infrare d flux from an aste roid (left, a radar 
shape model of Kleopatra by I Ostro et a/.L I2O0O1 and a “dirty im¬ 
age” (center) with atmospheric noise, obtained by transforming 
the incompletely sampled frequency plane (bottom, 6000 points, 
50 antennas). The smallest resolvable detail is approximately 
10 mas. Provided enough observations at different geometries are 
available (here eight observation runs were used), a shape model 
(right) can be constructed directly from the raw interferometric 
data. 


NIRVANA instrument can - in principle - be able to re¬ 
solve binary asteroids whose components are separated > 
20 mas. This will allow splitting many binary asteroids. 

Another important source of interferometric data in 
the near future is the Atacama Large Millimeter Array 
(ALMA). It will provide resolution ^ 5 mas at 0.3 mm and 
a dense mesh of baselines, thus ena bling “ i magin g” of hun¬ 
dreds of asteroids in the main belt (Busch, 2009h . A useful 
feature of the multimodal inversion is that the raw ALMA 
(and any other interferometric) data can be used directly 
as the original Eourier transform: there is no need to re¬ 
construct the image estimate. What is more, the recon¬ 
struction of the overall shape is insensitive to inaccuracies 
and uncertainties in the thermal model used because the 
model is determined mainly by the boundary of the pro¬ 
jection, the distributio n of brightness inside the boundary 
is much less important iviikinkoski and Kaasalainen , 2014 ; 


Viikinkoski et a/.L l2015h . An example of inversion of simu¬ 
lated ALMA data is shown in Eig. |9] 

Long-baseline interferometry can also represent a novel 
approach to determine the masses and the densities of aster¬ 
oids in a range of sizes and distances never studied before. 
In particular, modern interferometers such as the VLTI, 
LBTI and the Magdalena Ridge Observatory Interferometer 
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can spatially resolve binary asteroi ds discovered by p hoto- 
metric lightcurves in the main belt iCarrv et aZ.Ll2015h . The 
separation of the components of these systems is too narrow 
for traditional observational techniques such as AO at 10-m 
class telescopes. 

3.6. Direct size measurement with ESA Gaia 

The ESA Gaia mission, whose operations started in 
2014, will provide accurate astrometry and photometry of 
asteroids. Photometric data can be used the same way 
as sparse photometry from ground-based observatories 
(Sect. 13.11). The poten t ial of Gaia-like data was demon¬ 
strated by Cellino et al\ ( 20091) on data from the Hipparcos 
satellite. Apart from disk-integrated brightness, Gaia will 
provide also direct measurements of asteroid sizes, in a way 
much similar to the interferometry. Being designed as an 
astrometry mission at the /ras level, the PSF of Gaia is ac¬ 
curately known and stable. Deviations of the PSF from that 
of a point-like source can thus be measured and used to 
estimate the apparent size of asteroids. 

Owing to the amount of sources observed by Gaia, the 
satellite does not download images to Earth, but only 1- 
D flux profiles, corresponding to small 2-D windows cen¬ 
tered on targets, stacked along one direction, similar to in¬ 
terferometry in that respect. Interpretation of this 1-D pro¬ 
file relies therefore on the a priori knowledge of the 2-D 
flux distribution on the plane of the sky. In the creation of 
Gaia catalog, this will be done iteratively, starting from sim¬ 
ple spheres at zero phase angle to Anally use the spin and 
tri-axial elli psoid solutions determined otherwise from Gaia 
photometry ( Misnard et al\. 2007 ; 


Cellino et al. 


ume; Delbo and Harri]i 2002t Harris and Laserro^ 2002). 
Ground based telescopes can only observe in specific win¬ 
dows of the electromagnetic spectrum where the atmo¬ 
sphere is relatively transparent: i.e., the L, M, N, and Q 
bands at 3, 5, 10, and 20 pm, respectively. However, such 
observations are strongly affected by the variability of the 
transparency of the atmosphere and its thermal background. 
The background also receives contribution from the tele¬ 


scope and the optics. Hence, the thermal infrared obser¬ 
vation of asteroids from the ground is limited to relatively 
bright asteroids (V < 18 mag with 10-m class telescopes). 
Absolute calibration of the flux is rarely better than 5-10%. 
On the other hand, in space the instrument calibration is 
usually stable and there is no need to reduce the thermal 
background from the atmosphere. As a consequence, space 
based telescopes such as Spitzer can observe much fainter 


and smaller asteroids iMommert et al\. l2014h . with uncer¬ 


tainties in the calibration that can reach ^1% error. Also, 
from space the range of the observational wavelength is 
limited only by the detector technology, typically ^3.5- 
50 pm. At longer wavelengths, the telescope optics require 
cooling and the observation of faint objects is confronted 
with the background from solar system dust cloud and in¬ 
frared cirrus. 

In the Asteroids III era, the main source of thermal in¬ 
frared observations of asteroids has been the IRAS Mi¬ 
nor Planet Survey that collected obser vations of more 
than 2200 asteroids ( Tedesco et al\. 2002h . Since then, the 
NASA Wide-fleld Infrared Survey Explorer (WISE) has 
observed more tha n 130,000 main-belt asteroi ds (Mainzer 


2012ll201lh. a bout 500 


et al, this volume: \Masiero et ali _ _ 

near-Earth asteroids ( Main^^r^eta^ 20n, 2012bll3) . about 
1100 Hilda asteroids ( Gray et aTT 2012a), and almo st 2000 


Jupiter-Trojan asteroids ( Gray et al. . 2012bl 2011 ) in four 
infrared wavelengths at 3.4,4.6,12, and 22 pm-, the AKARI 
space telescope observed more than 5000 asteroids dur¬ 


ing it s mission UJsui et alx. 1201 iL 120131 : \Haseeawa et al 


2013h : the Spit zer space telescope observed hund r eds of 


2007 ). asteroids (e.g Emery et al . 20061: Trillins et ali 2010l : 


These measurements represent a great opportunity for 
multi-data inversion algorithms: the measurement of the ex¬ 
tension is direct, but clearly requires a realistic description 
of the projected shape on the plane of the sk y. Considering 
Gaia s pecifications and observing geometry. iMignanf et all 
( 2007h have estimated the fraction of asteroids for which 1- 
D dimension will be measured with a precision better than 
10%. This fraction is highly dependent on diameter, and 
20% of asteroids between 20-30 km will be measured at 
least once, while all asteroids larger than ^80km will be 
measured repeatedly. 

3.7. Disk-integrated radiometry 

Measurements of asteroids emission in the thermal in¬ 
frared - in general at wavelengths between 4-5 and 20- 
30 pm - are mostly used to determine the sizes of these 
bodi es {Mainzer et al, this volume; Delbo et al, this voT 


Licandro et alx. 1201 21 among others); the Herschel Space 


Observatory, that, due to its longer wavelengths, spanning 
55-671 iim, was primarily u sed to observe trans-Neptunian 


objects { Muller et al\. l201^ . For a review about all these 


missions and their results see Mainzer et al., this volume. 

The thermal infrared spectrum of asteroids carries in¬ 
formation about their size and surface properties, such 
as the thermal inertia, roughness and emissivity. These 
properties are typically derived by interpreting thermal in¬ 
frared data by means of thermal mode ls (Mainzer et al. , 
this volume; Delbo et al, t h is vol ume; {Delbo and Harrisi 
2002; Harris and Laserr^ 2002h . The “simple” ther¬ 


mal models that assume a spherical shape, a Lambertian 
emission of the surface, and a simplified calculation of 
the surface temperature distribution are used when we 
lack knowledge of the asteroid global shape, spin vector 
and rotation period, which is the majority of the cases. 


Thermal Model (NEATM. Harris 

. 1998; Delbo and Harris. 

2002: Harris and LaserroA\200il 

), or the Standard Thermal 

Model ISTM. {Lebofskv et ali 1 19861: Harris andLaserros. 

2002: Delbo and Harris 

o 

o 

CN 

c 


However, in order to derive the thermal inertia of an as¬ 
teroid from measurements of its thermal infrared emission, 
more sophisticated models, called thermophysical models 
(TPMs), are needed ( Spencer . 19901 : Spencer et al. [ 119891: 
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calculate the temperature distribution over the body’s sur¬ 
face as a function of different parameters, including the 
thermal inertia. In these models, the asteroid shape is usu¬ 
ally fixed and is modeled as a mesh of planar facets. The 
temperature of each facet is determined by numerically 
solving the one-dimensional heat diffusion equation using 
assumed values of the thermal inertia, with the boundary 
condition given by the time-dependent solar energy ab¬ 
sorbed at the surface of the facet (see Delbo et al, this 
volume, for a review). This latter quantity is calculated 
from the heliocentric distance of the asteroid, the value as¬ 
sumed for the albedo, and the solar incident angle. Macro¬ 
scopic surface roughness is usually modelled by adding 
hemispherical section craters of variable opening angle and 
variable surface density to each facet. Shadowing and mul¬ 
tiple reflections of incident solar and thermally emitted ra- 
di ation inside crat e rs are taken i nto account as described 
by Spencer ( 1990l): Emery et al. ( 1998 ): Rozitis and Green 
(1201 lb and xLaeerroa (119981). Heat conduction i s also ac¬ 
counted for wi t hin craters iSoencer et al\. 19891: Spencer , 
Il990l: Laeerro^ 19961: Delbo. 20041) . Surface roughness can 
be adjusted by changing the opening angle of the craters, the 
density of the c rater distribution, or a combination of both 
{ Mueller . 2007h . The total observable thermal emission is 
calculated by summing the contributions from each facet 
visible to the observer. Model parameters are adjusted until 
the best agreement with observational data is obtained, i.e., 
the least-squares residual of the fit is minimized, thereby 
constraining the physical properties (albedo, size, macro¬ 
scopic roughness, and thermal inertia) of the asteroid. 

From the point of view of multi-data inversion how¬ 
ever, the optimization of thermophysical parameters as de¬ 
scribed above is a two-step process - first, the spin and 
shape model is derived from one data type (photometry, 
radar,...), then this model is fixed and used for deriving 
thermophysical parameters from another data type (thermal 
infrared). This approach lacks the possibility to weight the 
two data types with respect to each other. Moreover, the 
thermophysical parameters can be very sensitive to small 
modifications of the input shape and spin, so various modi¬ 
fications of the shape should be tested to see how stable the 


solution is iHanus et al.i l2015h . Ideally, one should model 


shape and spin parameters together with thermal parame¬ 
ters. This multi-data approach using lightcurves and ther- 
mal infrared data simultaneously was succ essfully tested by 


Durech et al. ( 2012h : Durech et al\ ( 20141) and in principle 


it can be used also to data that are sparse in time. 


3.8. Radar 

Radar obseryations that measure the distribution of echo 
power in time delay and Doppler frequency (so-called 
range-Doppler or delay-Dopp l er me asurements) are dis¬ 
cussed in detail in \Ostro et al\ (l2002h and in Benner et al. 
(this yolume). The delay-Doppler projection is many-to- 



Fig. 10.— Example range-Doppler radar im ages (left) of asteroi d 
2000 ET 70 observed at Arecibo observatory iNaidu et al. I I20T3I) . 
Range and frequency resolutions are 15 m and 0.075 Hz, respec¬ 
tively. Range increases towards the bottom, Doppler frequency 
increases to the right. The simulated images (right) correspond to 
the reconstructed shape model in Fig. 


one mapping of a 3-D surface of the target into a 2-D “im¬ 
age”. Each pixel on the image represents a bin containing 
integrated echo power from surface elements that have the 
same distance from the radar and the same relative speed 
(due to the rotation of the asteroid). From the point of view 
of inversion, images in the range-Doppler plane are gener¬ 
alized projections that can automatically be handled with 
the general procedure discus s ed in Sect. |3and in detail in 


Viikinkoski and Kaasalainem (120141) and Wiikinkoski et al. 


120151) . In this approach, the multi-mode reconstruction is 
tuned to produce models with intermediate scale resolution 
(^ 1/10 of the diameter) since these are computationally 
inexpensive (can be obtained in a few minutes with a lap¬ 
top), and data sources other than radar do not contain more 
detailed information. An example of radar range-Doppler 
data of asteroid 2000 ET 70 is shown in Fig. [10] and the 
corresponding reconstructed shape model in Fig. [TT] If de¬ 
tailed radar data are available, s uch model can then be fur¬ 
ther refined ( Naidu et al. . 2013h with the radar techniques 
described by Benner et al. (this volume). 

Due to the steep decrease of echo power with the dis¬ 
tance to the object (the fourth power of the distance), as¬ 
teroids achievable by current radar facilities Arecibo and 
Goldstone are only close-approaching near-Earth asteroids 
or the largest members of the main belt. A number of exam- 
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Fig. 11.— Shape model of the asteroid 2000 ET 70 reconstructed 
from Arecibo and Goldstone delay-Doppler radar images with 
ADAM (Sect. [HD. 


pies of what can be reconstructed from radar data is given 
by Benner et al. (this volume). 

3.9. Asteroid interior 

We briefly discuss here the interesting possibility of 
peer into an asteroid. This is somewhat separate from 
the remote-sensing framework as the data can only be ob¬ 
tained in situ; on the other hand, such data can be acquired 
with the future space missions. The most practical op¬ 
tion is to plant radio transmitters/receivers on the surface 
of a kilometer-sized asteroid and measure low-frequency 
(~ 100 MHz) signals between these and an orbiter when 
they pass through the interior of the target. 

The m ost robust observables are simply s ignal travel¬ 
time data i Pursiainen and Kaasalainen . 2013h . These al¬ 
low an efficient formulation of the inverse problem via 
the refraction index and are relatively insensitive to noise 
and model error. Nevertheless, they suffice to give a 
coarse-scale picture of the general distribution of permit¬ 
tivity inside the asteroid, as well as the locations and 
sizes of large anomalies (sudden low- or high-density re¬ 
gions such as voids or heavier minerals). This approach 
has also been robustly tested in l aboratory conditions 
( Pursiainen and Kaasalainen . 2014bh . A more refined pos¬ 
sibility is to measure changes in the pulse profile, although 
this is more prone to errors ( Pursiainen and Kaasalain^ 
2014ah . The interior of the asteroid is practically impossible 
to model accurately in three dimension since it is supposed 
to have a number of cracks, voids, discontinuities etc., all 
refracting and reflecting the radio waves in complicated 
ways. Thus a very robust scheme is essential for extracting 
the available information with stability. Regardless of the 
data type, one or two transmitters on the surface are not 
sufficient for a unique solution. A tetrahedral configuration 
of four transmitters would be ideal, but this places heavy 
demands on the payload design. 


3.10. Extension of the model 

So far, we assumed that the asteroid can be described as 
a solid single body with constant spin vector, i.e., rotating 
along the axis with the maximum moment of inertia with 
a constant rotation rate. Although this model represents a 
typical asteroid, there are other configurations that can be 
also treated with an extension of the simple model. 
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Fig. 12.— Shape model of (90) Antiope seen form its equator 
(bottom) reconstructed from lightcurves and two occultations ob¬ 
served in 2008 and 2011 (top and middle). The model silhouette 
is projected on the plane of the sky, the solid lines are positive 
chords, the dashed ones are negative observations. 


3.10.1. Binaries 


Binary asteroids for m a significant part of the popula¬ 
tion of small asteroids. Pravec et al. ( 2006h estimated that 
the fraction of binaries is 15% for the near-Earth population 
and a similar fraction is assumed for the main-belt popula¬ 
tion in the same size range (Margot et al, this volume). The 
formation, dynamics, and physical properties of binary and 
multiple asteroid systems are discussed in detail in other 
chapters of this volume (Walsh et al., Margot et al). In gen¬ 
eral, modeling of such systems is more complex because 
the dynamics has to be taken into account. However, in 
some cases described below, the technique used for model¬ 
ing single asteroids can be used also for multiple systems or 
at least their primary components. 

For binary (or multiple) systems where the primary com- 
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ponent is much larger than the satellite, the photometric 
signal from the satellite can be neglected and the primary 
can be modeled as a single body. For such systems, the 
shape of the primary and the constraints on its gravitational 
quadrupole J 2 from the orbit analysis of the secondary 
can be used to investigate the distribution of the d ensity 


served during the occultation in 2011 \Colas et al\ . l2012h 
portray the two components to details unattainable by 
any other observational technique and the large set of 
lightcurves observed over many apparitions constrains the 
rotational parameters. A similar model can be obtained also 
by using lightcurves separately to create a scale-free model 
that is then scaled by occultation data (t he Shaping As¬ 
teroid s with Genetic Evolution algorithm, \Bartc7ak et al. , 
20141) . However, this two-step approach lacks the advan¬ 


tages of simultaneous inversion where the two data types 
can be weighted with respect to each other. 

3.10.2. YORP effect 

As described in detail by Vokrouhlicky et al. (this vol¬ 
ume), rotation state of small asteroids is affected by the 
anisotropic recoil of scattered sunlight and thermal radi¬ 
ation, which causes a net torque called the Yarkovsky- 
O’Keefe-Radzievskii-Paddack (YORP) effect. This effect 
secularly changes the obliquity of the spin vector and the 


rotation period. Whereas the former is too small to be mea¬ 
sured with current data, the latter has been measured on 
several asteroids (see Table 3 in Vokrouhlicky et al., this 
volume). If the change of the rotation period is larger than 
the uncertainty of the period, the change can be traced from 
apparition to appari tion as was the case for (54509) YORP 


iBerthier et al.l\20l4i\Takahashi and 5c/ieere.dl2014l) . 

For systems with comparable sizes, the problem be¬ 
comes complicated when the system is asynchronous, i.e., 
when the rotation period of the primary is different from the 
orbital period of the seco ndary. Such systems are usu ally 
modeled as two ellipsoids iScheirich and Praveo. 120091) . al¬ 
though more general models were cr eated from radar ob ser- 
yation (1999 KW 4 , for example, see l(9.yfro et all 2006 ). 

Fully synchronous binaries can be approximated by sin¬ 
gle bodies if the separation of their components is not 
large. Eyen a conyex model can proyide a good fit to the 
lightcuryes i Purech and KaasalaineA 2003h . Such model 
does not represent the true configuration of the system, of 
course, but it proyides the correct rotation period and direc¬ 
tion of the orbital plane. 

When the separation of components is larger, the system 
has to be modeled as two-component. Howeyer, from the 
modeling point of yiew, it is just a moderate modification 
of the nonconyex problem, where the system is described 
by only one rotation/orbital period and orientation of the 
normal of the orbital plane (parallel to the spins of the bod¬ 
ies). If the model is based on lightcuryes only, the spin 
and period parameters can be reconstructed accurately, but 
the uncertainty in sh apes is large. As has been shown by 
Marchis et al\ ( 201^ on (624) Hektor, the distinction be¬ 
tween a highly nonconyex single body, two bodies in con¬ 
tact, or two bodies orbiting each other is difficult to make. 

An example of reconstruction of a doubly synchronous 
binary system (90) Antiope from lightcuryes and occulta- 
tions is shown in Fig. [121 where the model is shown together 
with the silhouettes from occultations. Tens of chords ob¬ 


{ Lowry et al. . 200 ?!) . In other cases, the effect was much 


smaller and it reyealed itself by the discrepancy between 
the d ata and the model assuming the period to be con¬ 


stant iKa asalain en et al. 

Lowry et al . 20141) . 


2007 : durech et al\. 2008 . 2012 : 


The YORP effect is easy to include into the model. We 
assume that the rotation rate w changes linearly in time t 
as dw/df = V. Then the parameter v is another free pa¬ 
rameter of the modeling. Because the shift in the rotation 
phase increases quadratically in time as (j> = uit + l/2vt'^, 
eyen small changes do; /dt of the order 10 “® rad/d“^ can be 
detected with data sets coyering tens of years. 

In principle, the measured yalue of v can be compared 
with the theoretical yalue computed from the spin state, 
shape, size, and the thermal parameters of the surface 
with the density as a free parameter. Howeyer, due to 
sensitiyity of YORP on small scale details of th e shape 
that are far below the resolution of th e model LStatler . 
20091 : Kaasalainen and Nortunen ^ 2013 ), and the p r oblem 
of transyerse heat diffusion { Golu bov and Krusb\ . 2012 : 


Golubov et al\. 2014; Sevecek et al . . 2015 ). this can hardly 


be more than a rough comparison. 

Because the YORP effect scales as inyerse of the square 
of the size of asteroid, it becomes more important for small 
bodies, where it might be necessary to include it into the 
modeling if the data coyer a wider span of time. YORP is 
assumed to play an important role in many dynamical pro¬ 
cesses - the distribution of rotation periods and spin obliq¬ 
uities o f small asteroids i Pravec et al . 2008t Hanus etal\. 
2013bl) or the creatio n of asteroid binaries and asteroid pairs 
\Pravec et al . , 2010l) . for example. It is important to haye 
more asteroids with YORP detection. Better statistics of 
yalues of period change will help to constrain theories of 
YORP eyolution of small asteroids. 

3.10.3. Excited rotation 

Asteroids rotating in the relaxed mode are fully de¬ 
scribed by the spin axis direction, rotation rate, and the 
initial orientation. Howeyer, some asteroids are in an ex¬ 
cited rotation state, which can be described as a rotating free 
top. Mo re parameters are ne eded to describe this tumbling 
motion iKaasalainen 120011) . The reason why some aster¬ 
oids are in this state can be (i) primo rdial, (ii) collisional 
excitation iHenych and Pravea 1201 3h . or (iii) end state of 
YORP-driyen spin down. The approach to the modeling is 
in principle the same as for asteroids in principal axis rota¬ 
tion, only the orientation for a giyen time is giyen by solying 
differential eq uations. From lightcur yes, models of aster- 
oids 2008 TC^ iScheirich et ~^. , 2010l) and (99942) Apophis 
iPravec et all 120141) haye been deriyed. A model of aster- 


13 




















































































































90 - 


30 - 


- 30 - 


- 90 - 




• •• 


10 


Size [km] 


100 


Fig. 13.— Distribution of pole ecliptic latitude /3 with respect to 
the size of the asteroid for ~ 800 asteroids. The the scale on the 
vertical axis is linear in sin /3, which makes the vertical distribu¬ 
tion of points in this plot uniform for an isotropic distribution of 
spin vectors. 


ical properties can be drawn. The second approach is to 
concentrate on selected targets, obtain as many different 
data types as possible, create detailed models of these aster¬ 
oids, and extrapolate the obtained results to the whole class 
of similar objects. We describe in this section new research 
areas that directly benefit from availability of spin solutions 
and 3-D shape models. 

4.1. Spin-axis distribution and evolution 

One of the main results of the lightcurve inversion is 
the increasing sample of asteroids with known orienta¬ 
tion of the spin axis. For main belt asteroids, the long- 
kno wn lack of asteroids with poles close to the eclip¬ 
tic ( Kr\s7C7\nska et all 2007 ) was confirmed and it was 
shown that i t is m ore pronounced for smaller asteroids 
Wanus et ali 1201 ih . In Fig. [13] we plot the distribution 
of pole latitudes for ^ 800 asteroid models with respect 
to the ir size (an updated version of Fig. 5 in Hanus etal\. 


oid (4179) Toutatis iHudson et a7ll2003l) was derived from 
radar data and lightcurves. This model was later compared 
with the fly-by images of Chang’E-2 mission - the general 
shape was in agreement , altho ugh there were some minor 
discrepancies \Zou et a/.Ll2014l) . 

3.11. Procedures: Convexinv, KOALA, and ADAM 

Software for the inverse problems is available at DAMIT 
(http://astm.troja.mjf.cuni.cz/projects/asteroids3D). Convex 
is a procedure for lightcurve inversion, while ADAM (All- 
Data Asteroid Modeling) is a collection of functions from 
which o ne can tailor an i n versio n procedure for any data 
sources iViikinkoski et al.i l2015h . An earlier version of 


201 Ih . The size-dependent structure can be explained by 


this is called KOALA (Knitted Occultation, Adapti ve op¬ 
tics and Lightcurve Analysis, \Carrv et all l2012h : this 
is based on lightcurves and silhoue tte contours obtain¬ 
able from images and occultation s { KaasalaineA l201 it 
Kaasalainen and ViikinkoskA\20l'2xi . KOALA is especially 
suitable for lightcurve and occultation data, while ADAM 
allows the use of any images (camera, radar, or interfer¬ 
ometry) with or without lightcurves without having to pro¬ 
cess them to extract contours or other information. Both 
KOALA and ADAM can be used for lightcurves only but, 
as discussed earlier, this is not reliable. Any resulting 
shape should only be taken in the global sense (as a more 
realistic-looking rendering of a convex solution), and the 
details or nonconvex fe atures are seldom likely to be real 
( Viikinkoski et al. . 2015h . 


4. WHAT HAVE WE LEARNED FROM 3-D SHAPES? 


the YORP effect that is more effective on smaller asteroids 
(<30km) and pushes them into extreme values of obliq¬ 
uity (0° or 180°). This corresponds to the clustering of pole 
latitudes towards values of ±90°. Although there are ob¬ 
servation and modeling biases that affect the distribution 
of poles in the sample of available models, their effect is 
only ma rginal compared to the strong anisotropy seen in 
Fig. [T3] ( Hanus et all 2011 ). The spin-axis orientation of 
W6h smaller asteroids (< 5 km) is still not known due to the 
lack of models. For the largest asteroids (> 60 km), there is 
a statistically significant increase of prograde rotators (98 
prograde v s. 63 retrograde in Fig. T^, pr obably of primor¬ 
dial origin { Johansen andLacerda . 2010h . 

A different appro ach to the problem o f spin-axis distri¬ 
bution was used by Bowell et al\ ( 2014 ). They analyzed 
variations of the mean brightness with the ecliptic longi¬ 
tude, from which they estimated ecliptic longitudes of spin 
axis for about 350,000 asteroids and revealed a clearly non- 
uniform distribution. However, the explanation of the cause 
for this non-uniformity is still missing. 

With increasing sample of models, it is also possible to 
study the distributio n of spin axes of mem bers of collisional 
families. Results of Hanus etai\ ( 2013a ) agree with theo¬ 
retical expectations; if the spread in proper semimajor axis 
increases with decreasing size due to the Yarkovsky effect 
(Nesvorny at al., Vokrouhlicky et al, this volume), asteroids 
closer to Sun than the center of the family should rotate ret¬ 
rograde, those farther should rotate prograde. This behavior 
is shown in Fig.[T4|for Flora family. 


The purpose of modeling methods described in previous 
sections is to reveal new facts about the nature of asteroids. 
The approach is made on two fronts: The first strategy is 
to use the most abundant data sources (photometry in vi¬ 
sual and thermal IR) to produce many low-resolution mod¬ 
els that will be a statistically significant sample of the whole 
asteroid population from which conclusions about the phys¬ 


4.2. A pre-requisite in many cases 

Unlike most astronomical objects, the viewing geometry 
of asteroids continuously changes due to their motion rela¬ 
tive to observer. As a corollary, knowing their rotation pe¬ 
riod and spin-vector coordinates is crucial to interpret pro¬ 
jected size measurements, and to tie together observations. 

As already described in Sect.|2 the spin and 3-D shape 


14 















































Fig. 14.— Distribution of pole ecliptic latitude (5 with respect to 
the proper semimajor axis Up of the Flora family members. The 
relative size of each circle corresponds to the size of the asteroid. 


are required to interpret thermal infrared radiometry, stellar 
occultations, or apparent sizes measured by Gaia without 
biases. Similarly, there have been long running discussions 
to explain inconsistent spectral measurements of asteroids, 
that can be easily solved once the spin properties (period 
and orientation) are k nown: see, for instance, the discus¬ 
sion o n (83 2) Karin by Sasaki et al. (2004), Vemazza et al\ 
( 2007h. and Chapman et al ( 2007 ). or on (21) Lutetia by 
Barucci et al. ( 20121) . 


4.3. Density, composition, and internal structure 

Density is one of the most fundamental properties to 
constrain the composition of asteroids and investigate their 
internal structure {Scheeres et al, this volume). With the ex¬ 
ception of binary asteroids with observable mutual eclipses 
{Margot et al, this volume), both mass and volume are re¬ 
quired to determine the density of an asteroid. Estimating 
the mass of an asteroid by measuring its gravitational influ¬ 
ence on other objects is a challenge because of the relative 
low mass of asteroids compared to other planetary objects. 
The number of mass determinations thus limits the number 
of density estimates. Although there are diameter, hence 
volume, estimates for all asteroids with a mass determina¬ 
tion, the uncertainty in volu me generally d ominates the bal¬ 
ance on density uncertainty (ICarnil2012l) . 

In this respect, 3-D shape models are required to deter¬ 
mine accurate volume. The level of potential biases in¬ 
creases with stronger assumptions on the shape, and ac¬ 
curacy accounting for systematics improves from sphere, 
to ellipsoid, to convex hull, to the real shape. In the 
decade since Asteroids III, the number o f density es t imates 


has increase d from 20 to 300 objects iBritt et al\. 12002 : 


ICarni 20I2h . Among these, the most reliable are derived 
from binary systems for which the volume of the primary 


2006t Sheoard et al , 

2006; 

Descamos et al . 2009 

, 2011 

Fan 2 et al . 2011; Marchis et al , 20131). 


A detailed description on the density of asteroids and 


their internal structure can be found in Scheeres et al. (this 
volume), from under-dense asteroids, hosting large voids, 
to over-dense asteroids, likely differentiated. In the con¬ 
text of this chapter, it is important to highlight than only 
density estimates more precise than 10-20% can be used to 
discriminate between different analogue meteorites and can 
provide insights on the internal structure. Such accuracy 
can only be achieved with volume known to 5-10% or bet¬ 
ter, which means that a proper description of the 3-D shape 
is needed. 


4.4. Determination of surface properties by means of 
TPMs 

The spin state and shape model of an asteroid is in¬ 
put information for the TPMs. While in the Asteroids III 
era, shape and spin information were available only for a 
handful of asteroids, preventing application of TPMs to a 
large number of these bodies, this situation has drastically 
changed in the last few years. Physical properties (such as 
the value of the surface thermal inertia) of about 60 aster¬ 
oids are now available thanks to the application of TPMs 
(see Delbo et al, this volume). A remarkable improve¬ 
ment in this held is also represented by the availability of 
high quality thermal infrared data as those produced but the 
WISE and the Spitzer space telescopes. In the next few 
years we expect the number of asteroids with known ther¬ 
mal inertia values to grow thanks to the availability of more 
shape and spin state models. 


4.5. Surface re-arrangement 

Erom the spin, 3-D shape, and density, the local gravity 
at the surface can be computed. Unsuspected physics has 
been unveiled with the modeling o f the ne ar-Earth asteroid 
(66391) 1999 KW4 bv lO^tro et al\ (l2006h . Some small as¬ 
teroids present an equatorial bulge, presumably generated 
by regolith migration toward lower gravity regions. This 
process can even form binary systems if the asteroids spin 


fast enough malsh et a/.Ll2008t I j/arm et a/.L 120091) . 


4.6. Cratering events 

An evident outcome of shape modeling is the capa¬ 
bility to detect large impact craters and basins. Aside 
from the spacecraft encounters, the first detection was the 
large impact basin o n Vesta, progenitor of the Vestoids 
( Binzel and xA 1993 ), detected with the HST ( Thomas et al\. 
2997 ) and confirmed by NASA Dawn spacecraft Wussell et al 
2012h . Another case is the recent impact suffered by aster- 
oids (596) Sche i la, de tected by the presence of a dust tail 


iBodewits et a/.L 1201 Ih . Lightcurves obtained before and 


after the impact, under similar geo metries are different, re - 
vealing different surface properties iBodewits et g/.L 120141) . 
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Table 1 

A LIST OF OBSERVATION TECHNIQUES AND DERIVABLE PHYSICAL PROPERTIES. 


Technique 

Period 

Spin 

Size 

Shape 

Thermal Inertia 

Number of models 

Asteroids III Asteroids IV Asteroids V 

Photometry 

X 

X 


X 


30 

500 

10"* 

Images 


X 

X 

X 


5 

50 

10^ 

Occultation 


X 

X 

X 


5* 

50 

10^ 

Radar 

X 

X 

X 

X 


10 

30 

10^ 

Radiometry 



X 


X 

10* 

20 

10-* 

Interferometry 


X 

X 

X 


5t 

< 10 

10^ 

Fly-by 

X 

X 

X 

X 

X 

6 

10 

< 15 


The “X” mark indicates which physical properties are derivable from which technique. The number of models 
available at the time of Asteroids HI book (Asteroids III) and now (Asteroids IV) is only approximate. The Asteroids V 
column is an order-of-magnitude estimate for the next decade. 

‘Ellipsoidal models. 

tHST/FGS. 


4.7. Mass distribution 


Because some asteroids are less or more dense than their 
most-likely constituents, the question of the mass distribu¬ 
tion (denser material or voids) in their interior can be asked. 
This question is intrinsically tied with the study of the grav¬ 
ity field around the asteroid. Th e latter has been me asured 
during spacecraft encounters (see Miller et ^ 2002 , for in¬ 
stance), but studies from Earth-bound observations have re¬ 
cently appeared. By comparing the spherical harmonics of 
the gravity field as determined from the orbit of a natural 
satellite, with the expected coefficients resulting from the 
3-D shape model, the h ypothesis of homogeneous mass dis 
tribut i on can be tested \Vac hier et al. 20121; Berthier et al. 


1201 4 \Takahashi and Scheere^YlDX^i ). 


5. FUTURE 


In the decade since Asteroids III, where the principles 
of lightcur ve inversion based on dens e-in-time series were 
presented ( Kaasalainen et ai\ 2002bh . the number of mod¬ 
els has seen a tenfold increase - from a few tens to a few 
hundreds (Table [T]i. The increasing availability of sparse- 
in-time photometry, and its appropriate handling in the in¬ 
version, coupled with a dramatic increase of computer time 
(thanks to projects like Asteroids®Home) have made this 
possible. Upcoming all-sky surveys such as Pan-STARRS, 
LSST, and Gaia are expected to produce enormous data 
sets and there is little doubt that thousands of models will 
be derived in the next decade. Our knowledge on non- 
gravitational effects such as YORP and Yarkovsky will di¬ 
rectly benefit from this larger sample. 

With the large number of data and models, new chal¬ 
lenges will arise - how to extract scientifically interesting 
information from a large set of models of asteroids? With 


big data flows, the processing has to be automated, with ef¬ 
fective data processing. The obvious search for correlations 
has to be done with care because of large biases in the set of 
models. Although the importance of detailed models of in¬ 
dividual well-studied asteroids will be important, the main 
shift in paradigm and probably the main source of interest¬ 
ing findings will be in tens of thousands of asteroid mod¬ 
els derived from photometry in optical and thermal infrared 
wavelengths. Only a few years ago, only photometry in the 
visible was available for a large number of asteroids. The 
situation drastically changed with WISE catalog of thermal 
fluxes for 150,000 asteroids. Automatized procedures ca¬ 
pable of dealing with photometry in the visible and ther¬ 
mal infrared will yield not only 3-D shape and spin state 
for thousands of asteroids, but also their diameter, albedo, 
and thermal inertia of their surface. The later being cru¬ 
cial in interpreting asteroid mineralogy once coupled with 
spectroscopy {Reddy et al, this volume). Understanding of 
observational and modeling biases will be crucial for cor¬ 
rect interpretation of the results. Connecting spin and shape 
distribution of asteroids with their orbital and spin evolution 
will hopefully lead to a clear picture of the evolution of the 
main asteroid belt. The ultimate goal here is the connec¬ 
tion of models of evolution of the main belt with spin-axis 
evolution and current distribution of asteroid physical prop¬ 
erties. 

Another approach that is complementary to modeling in¬ 
dividual objects, is modeling distribution functions of pa¬ 
rameters of interest. With any inversion technique and pho¬ 
tometric data quality, the number of models will be always 
much lower than the number of known asteroids just be¬ 
cause it takes time to collect enough data at different ge¬ 
ometries. So instead of aiming to create unique models 
for a statistically significant sample of the (sub)populations 
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with known observational and modeling bias, one can use 
other observables than the time-brightness pairs. With this 
approach, there is almost no “wasting” of data because es¬ 
sentially all photometric points are used in the statistical 
approach. The aim is to model characteristics of a given 
asteroid population as a whole when there are not enough 
data to model individual members. One of the observables 
can be, for examp l e, the mean brightness over one appari¬ 
tion i Bowell et al\. 2014 ) and its dispersion that correlates 
with lightcurve amplitude. If there are not enough data 
points to estimate mean brightness and its variance, the 
statistics of scatter of individua l pairs of observations can 
be used ( Szaho and /fii; J . 20081) . If the data come in pairs 
of observations separated by a constant time interval (like 
with Pan-STARRS), the rotation period can be estimated 
i Purech et ali 2007 ). 

Another challenge for next decade resides in the deriva¬ 
tion of a large sample of higher-resolution models, based on 
multiple data sources. Building inversion techniques capa¬ 
ble of handling different data sources had been highlighted 
in Asteroids III as the next step, and this became reality 
over the past few years. The models based on multiple data 
sources are more realistic, more reliable, more precise than 
the independent analysis of individual data sets. If practi¬ 
cal issues of merging data sets arise, mathematical solution 
have been exposed, and algorithms and software are ready 
and freely available. The current sample of models derived 
by such methods is nevertheless still limited. If the multi 
data analysis is more efficient than the independent analy¬ 
sis of the individual data sets, the practical problem of data 
sharing arises. Although the tools to invert multiple data 
sources are ready and freely available, the data are still hard 
to get and harder to share. To maximize the scientific out¬ 
put, it is necessary to combine all available data sources. 

The other observing techniques (i.e., disk-resolved 
imaging, stellar occultations, Gaia size measurements, or 
radar echoes) will never provide data for more than a few 
thousands asteroids. The apparent sizes by Gaia will how¬ 
ever be measured with a high precision, and will concern a 
decent sized sample. Because the data will be made pub¬ 
lic, they should be used on a regular basis for modeling. 
The number and quality of profiles derived from stellar 
occultations have always been restricted due to the intrin¬ 
sic complexity of prediction, and cost of equipment. The 
parallel availability of low-cost telescopes and cameras, 
together with the publication of Gaia stellar catalog and as¬ 
teroid orbits, will open a new age for stellar occultations. 
Any mass-production procedure of asteroid models should 
therefore be able to deal not only with photometry (visible 
and infrared) but also with stellar occultations. 

Finally, if the bulk of asteroids are single objects, rotat¬ 
ing along their shortest axis, a small fraction will either be 
in tumbling state, or have satellites. Automatic modeling 
procedures should be able to detect such cases, triggering 
detailed analysis of these peculiar targets. 

The rising interest of private companies for near-Earth 
objects as mining resources will inevitably lead to further 


development of analysis of remote-sensing data, building 
dedicated ground-based or space telescopes, and eventually 
moving from asteroids surface into their interiors. 
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